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Abstract 

We solve the linear advection-diffusion equation with a variable speed on a semi-infinite 
line. The variable speed is determined by an additional condition at the boundary, which 
models the dynamics of a contact line of a hydrodynamic flow at a 180° contact angle. We 
use Laplace transform in spatial coordinate and Green's function for the fourth-order diffusion 
equation to show local existence of solutions of the initial- value problem associated with the 
set of over-determining boundary conditions. We also analyze the explicit solution in the case 
of a constant speed (dropping the additional boundary condition). 

1 Introduction 

Contact lines are defined by the intersection of the rigid and free boundaries of the flow. Flows 
with the contact line at a 180° contact angle were discussed in [2l|l], where corresponding solutions 
of the Navier-Stokes equations were shown to have no physical meanings. Recently, a different 
approach based on the lubrication approximation and thin film equations was developed by 
Benilov &: Vynnycky [l]. 

As a particularly simple model for the flow shown on Figure [H the authors of [I] derived the 
linear advection-diffusion equation for the free boundary h{x,t) of the flow: 

dh d^h , ^dh , , 

The contact line is fixed at x = in the reference frame moving with the velocity —V{t) and is 
defined by the boundary conditions h\x=o = 1 and hx\x=o = 0. The flux conservation condition 
is expressed by the boundary condition hxxx\x=o = (take a'^ = 3 in equations (5.12)-(5.13) 
in IB). 

We assume that h,hx,hxx — )• as a; — t- oo: in fact, any constant value of h at infinity is 
allowed thanks to the invariance of the linear advection-diffusion equation (jl.ip with respect to 
the shift and scaling transformations. With three boundary conditions at x = and the decay 
conditions as x — )• oo, the initial- value problem for equation (jl.ip is over-determined and the 
third (over-determining) boundary condition at x = is used to find the dependence of V on t. 
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Figure 1: Schematic picture of the flow between rigid boundaries. 



We shah consider the initial-value problem with the initial data h\t=o = ho{x) for a suitable 
function Hq. In particular, we assume that the profile hQ(x) decays monotonically to zero as 
X — >• oo and that is a non-degenerate maximum of ho such that /io(0) = 1, /io(0) = 0, and 
/iq(0) < 0. If the solution h{x,t) losses monotonicity in x during the dynamical evolution, for 
instance, due to the value of hxx{0,t) crossing from the negative side, then we say that the 
flow becomes non-physical for further times and the model breaks. Simultaneously, this may 
mean that the velocity V{t) blows up, as it is defined for sufficiently strong solutions of the 
advection-diffusion equation by the contact equation: 

hxxxxx{0,t) = V{t)h^^{0,t), (1.2) 

which follows by differentiation of (jl.ip in x and setting x — )• 0. 

The main claim of [Ij based on numerical computations of the reduced equation (II. Ih as well 
as more complicated thin- film equations is that for any suitable ho, there is a finite positive time 
Iq such that V{t) — )• — oo and /ixz(0,t) — ?• —0 as 1 1 ^o- Moreover, it is claimed that V{t) behaves 
near the blowup time as the logarithmic function of t, e.g. 

y(t)~Cilog(to-t) + C2, as tt^o, (1.3) 

where Ci, C2 are positive constants. 

This paper is devoted to analytical studies of solutions of the advection-diffusion equation 
(jl.ip and the effects coming from the inhomogeneous boundary condition hxxx\x=o = ~^ associ- 
ated with the flux conservation. In particular, we rewrite the evolution equation for the variable 
tt = /la; in the form 

ut + u^xxx = V{t)ux, X > 0, i > 0, (1.4) 
subject to the boundary conditions at the contact line 

u\x=0 = 0, Uxx\x=0 = — Uxxx\x=0 = 0, t > 0, (1.5) 
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where the boundary conditions Uxxx\x=o = hxxxx\x=o = follows from the boundary conditions 
h\x=o = 1 and hx\x=o = as well as the original evolution system (jl.ip as j; — )• 0. 

To simplify the problem, we shall also consider the model for given constant V{t) = Vq and 
drop the third over-determining boundary conditions at the contact line: 

Ut + Uxxxx = VqUx, X > 0, t > 0, 

u\x=o = 0, t > 0, (1.6) 

'^xx\x=0 — 2' t ^ 0. 

Both problems (ll.4|) ~ ()1.5p and ()1.6p are considered under the initial condition u\t=o = uq{x) with 
no(0) = 0, «o(0) < 0, and Uo(0) = — ^, as well as the decay condition u,Ux,Uxx — )• as x — )• oo. 

Using Laplace transform in spatial coordinate and Green's function for the fourth-order dif- 
fusion equation, we derive an explicit solution of the boundary- value problem (|1.6p . In the case 
Vb = 0, we show that the inhomogeneous boundary condition hxxx\x=o = Uxx\x=o = — | leads to 
the secular growth of the boundary value hxx\x=o = Ux\x=o to positive infinity as t — )• oo. As a 
result, even if hxx\x=o < initially, the convexity of the solution h{x,t) at the boundary x = 
is lost in the finite time. In the case Vq < 0, we show that no secular growth is observed but 
the convexity of the solution at the boundary is still lost in the finite time. Applying the same 
method, we prove local existence of solutions of the original boundary- value problem (jl.4p - ()1.5p . 
This prepares us to tackle the original conjecture on the finite-time blow-up in the dynamical 
behavior of the model, which is still left opened for forthcoming studies. 

The remainder of this paper is organized as follows. Section 2 reports explicit solutions of 
the boundary-value problem (ll.6|) for Vq = and Vq ^ 0. Section 3 gives the local existence 
result for the original problem (|1.4p - (|1.5p . Appendix A lists properties of Green's function for 
the fourth-order diffusion equation. 



2 Solution for V{t) = Vq 

Because V{t) is nonconstant for the original problem (II. ip . the Laplace transform in time t is 
not a useful method for this problem. On the other hand, since the boundary-value problem is 
formulated for half-line, we can use Laplace transform in space x: 

/•oo 

U{p,t)= e-P''u{x,t)dx, p>0. (2.1) 
Jo 

We shall develop this method to solve the boundary-value problem (II. 6p . The explicit solution 
of this problem will help us to analyze the consequences of inhomogeneous boundary condition 
Uxx\x=o = and the constant advection term V{t) = Vq on the temporal dynamics of the 
advection-diffusion equation with the fourth-order diffusion. 
Let us denote the boundary values: 

/3(t) = Ux\x=0, l{t) = Uxxx\x=0- (2.2) 
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Using Laplace transform (j2.ip . we rewrite an evolution problem associated with the advection- 
diffusion equation (ll.6p : 

/ Ut+p^U-VopU = j{t)-lp + p^m, t>0, 

1 UU=o = Uoip), ^'-^^ 
where Uq is the Laplace transform of uq = u\t=Q. By using the variation of parameters, we obtain 

+ ^ ^-{t-s)p^+{t-s)v,p (^^^^ -\p + P^(3{s)^ ds. (2.4) 
Using the inverse Laplace transform in x, we write this solution in the form: 

-| pc+ioo / POO \ 

e-('-^^p'+('-'^^oP j^^(g) _ + p^m^ ds^ dp, (2.5) 



c+ioo 
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27rz 

where Re(c) > so that the singularities of the integrand in the complex p-plane remain to the 
left of the contour of integration. 

If t > is finite, uq e L^(R+), and /3,7 e Lg'^(M+), Fubini's Theorem implies that the 
integration in p and in y, s can be interchanged. Let us introduce Green's function Gt{x) for the 
fourth-order diffusion equation (see Appendix A): 

-1 pc-\-ioo -1 poo -1 poo 

^ I .VX-tp^^^ _ ±_ / ,~tk^+ikx.. _ i / -tk* 



Gt{x) = / ef^-*f dp = — e-"" +'^''dk = - / e""^ cos{kx)dk. 

2vri Jc+ioo 27r it Jq 

Using Green's function, we can rewrite the solution (j2.5p in the implicit form: 

noo 1 /■* , 

u{x,t) = / Gtix + Vot - y)uo{y)dy - - G't_,,{x + Vo{t - s))ds 
Jo ^ Jo 

+ [ [Gt-s{x + Vo{t - s))7(s) + G';_,{x + Vo{t - s))P{s)] ds. (2.6) 



The solution is said to be in the implicit form, because the functions f3{t) and 7(t) determined 
by the boundary conditions ()2.2p are not specified yet. 

We verify that lim^^^oo u{x, t) = 0, no matter what /3 and 7 are, as long as they are bounded 
function of t. Indeed, by the Lebesgue's Dominated Convergence Theorem, we have 



/•oo 

/ Gt{x + Vot — y)uo{y)dy ^ as x — )• 00 
^0 



if no G L^(R), because Gt{x) — )• as x — )• 00. On the other hand, the other three convolution 
integrals are bounded if /3, 7 G Lj^^(]R+) and t > is finite, because Gt, G[, and G'l have integrable 
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singularities at t = 0. By the same Lebesgue's Dominated Convergence Theorem, these three 
integrals decay to zero as a; — )• oo. 

The functions (3{t) and 7(t) are to be found from the integral equations obtained at the 
boundary conditions u{0,t) = and Ux{0,t) = f3{t). These derivations are performed separately 
for the cases of Vq = and Vq / 0- 

2.1 Case Vq = 

We rewrite the solution (|2.6p for Vq = 0: 



poo -I I't 

u{x,t)= Gt{x - y)uo{y)dy - - G[_,{x)ds + [Gt-s{x)i{s) + G'U{x)m] ds. (2.7) 

JQ ^ Jo Jo 

Using the boundary values ()A.3p and (|A.4p for the Greens function Gt{x) and the boundary 
condition u(0, t) = 0, we evaluate this expression at a; = and obtain an integral equation for j3 
and 7: 

To use the boundary condition Ux{^,t) = /3{t), we shall recall from equation (|A.5P that the 
function G'l'{x) behaves like 0{t~^) for any x > and hence is not integrable in f at t = 0. 
Therefore, we have to be careful to differentiate the solution in the above convolution form. The 
last term of the solution (|2.7p can be computed by using the Fourier transform: 



t 1 /"+00 / f* 



v{x,t) := J G'l_^{x)l3{s)ds = ^ J {ikYe"''' i^J e'^ ^'~')(3{s)ds ] dk. 
Differentiating this expression in x and integrating by parts in s, we obtain 

Vxix,t) = {ikfe"'''(^J^e-'''^'''^P{s)dsjdp 



^Tri ./-oo ^ \Jo ds 



-kHt-s)^ P^^S^ds] dp 



+00 ikx / i-t 



1 r'^ 



27rz ./-no k 
1 
2 



m-mHt{x)- 1 Ht-s{x)P\s)ds. (2.9) 



where 



^*(^) ^= ^ / = - / P^dk = / Gt{y)dy. (2.10) 

Here we note that all integrals are evaluated in the principal value sense, because the half-residue 
at A: = is canceled out in the resulting expression ()2.9p . Also we note that the decay of Vx to 
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zero as a; — 7- oo is satisfied because of the symmetry and normalization of Gt in ()A.6p . We can 
now use the boundary conditions (IA.4P and Ux{^-,t) = j5{t) to obtain the exact value for j3{t): 

m = 2/ G[{-y)u^{y)dy - / G^,(0)(is 
Jo Jo 

= 2 r G[{-y)u,{y)dy+^^^t^'\ (2.11) 
Jo ^ 

After I3{t) is found uniquely from ()2.1ip . ^{t) is found uniquely from the integral equation ()2.8h . 
This computation completes the construction of the exact solution of the boundary-value problem 
(II. 6p for Vo = 0. Now we turn to the analysis of obtained solution. 

Theorem 1 Consider the advection-diffusion equation il.6\} for Vq = with the initial data 
uq G L^(M+). Then, there exists a solution u G L°°(M_|- x M^) of the evolution problem in the 
explicit form ^2. 7| ), where /3, 7 E Lj^^(M+) are defined by h2. 8(j and h2.11\) with Ynat^oo P{t) = +00. 



Proof. The convolution integral in the explicit expression (|2.1ip can be analyzed from the 
Green's function (jA.Sp . If uq G i^M), then 



Gt{-y)uo{y)dy 



< 



\9 ||l°° IfoIIli 
iV2 



t > 0. 



Therefore, /3 G and /3{t) ~ t'^/^ as t ^ 00 due to the second term in (l2TT]l . Now, the 

integral equation (|2.8p for j{t) with a weakly singular kernel is well defined and solutions exist 
with 7 G L£;(R+). Similarly, the solution u G L°°(R+ x R+) is weh defined by ([231). □ 

Remark 1 One can show that there is no singularity of the solution for f3{t) as t ^ so that 
/5(0) = 1*0(0) by continuity. Also, one can show that the solution of the integral equation (E 
for ^{t) exists in the closed form: 7(i) = 2 Gt{—y)u'Q {y)dy . 

Coming back to the original question, if uo(0) = 0, tio(O) < 0, and Uo(0) 
a finite value of to G (0, 00) such that Ux\x=o > for all t > to, that is, h{x,t) loses monotonicity 
at X = in a finite time to (recall that u = h^)- This dynamical phenomenon occurs because of 
the inhomogeneous boundary conditions Uxx\x=o - 
in the fourth-order diffusion equation ()1.6p . 



^, then there is 



I even in the absence of the advection term 



2.2 Case Vq 

We have the solution in the implicit form (j2.6p and we need to derive integral equations on the 
unknown function /3(t) and 7(t). One integral equation follows again from the boundary condition 
n(0,t) = 0: 

- r [Gt^siVoit - s)Us) + G'UiVoit - s))/3(s)] ds 
Jo 

Gt{Vot-y)uo{y)dy-l [ G't^,{Vo{t - s))ds. (2.12) 

^ JO 
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To find another integral equation from the boundary condition Ux{0,t) = P{t), we have to use 
the technique explained in Section 2.1 and to compute the derivative of the solution ()2.6p in x: 

Ux{x, t) = / G'tix + Vot - y)uo{y)dy - - / G'l_,{x + ¥^{1 - s))ds 
Jo ^ Jo 

+ [ G't_,{x + Vo{t - s)Ms)ds + - m)Ht{x + V^t) 

Jo ^ 
t l-t 

Ht-s{x + Vo{t- s))(3'is)ds + Vo Gt-six + Vo{t-s))(3{s)ds. (2.13) 



m + 2mHt{Vot) + 2 r Ht-siVoit - s))(3'{s)ds 

Jo 
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We can now use the boundary condition Ux{0,t) = (3{t) to obtain another integral equation for 
(5 and 7: 

ft 

'0 

-2^0/ Gt-s{Vo{t-s))P{s)ds-2 f G[_,{Vo{t - s))^{s)ds 
Jo Jo 

f'OO rt 

= 2 G'tiV^t - y)uo{y)dy - G'UiVoit - s))ds. (2.14) 
Jo Jo 

The system of integral equations (|2.12p and (j2.14p completes the solution (j2.6p for the case Vq 7^ 0. 
Because of the original motivation to study behavior for large negative V{t) in ()1.3p . we shall 
analyze the obtained solution for Vq < 0. 

Theorem 2 Consider the advection-diffusion equation il.6\} for Vq < with the initial data 
uq G L^(M_|_). Then, there exists a solution u £ L°°(W^ x M_|.) of the evolution problem in the 
explicit form \2. 6\) . where /3,7 G L°°(]R_|.) are defined by 112.1^) and \2. 14\ ) with 

Proof. Similarly to the proof of Theorem [H it is easy to show from the integral equations 
(IXT^ and ([XTip that if uq £ L^{^+), then /3,7 G L^^{R+). We shall now compute the limit of 
/3(t) and 7(t) as i — )■ 00: 

/3oo := lim /3(t), 700 := lim j{t). (2.16) 

To deal with the first integral equation ()2.12p . we first notice the explicit computation by using 
the Fourier transform: 

fit) := f G[^,mt-s))ds 
Jo 



fOO / ft 

27r 



2tt ./__ - iVr 
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where the integrals in s and k can be interchanged by Pubini's Theorem and the integration is 
performed in the principal value sense. We can now explicitly compute the limit as f — )• oo by 
using Lebesgue's Dominated Convergence Theorem: 

hm /(t) = — / dk 



t^oc'^ ^ 2tt A;3 - iVo TT Jo k^ + 3\Vo\^/^ ' 

This computation gives the last term of the integral equation ()2.12p as t — )• oo. To deal with the 
first term on the right-hand side of (|2.12p . we write 

/"OO 1 roo / foo \ 



oo 



where 



By Lebesgue's Dominated Convergence Theorem, this integral converges to zero as t — )• oo as 
long as uo e L^(M+). 

To deal with the second term on the left-hand side of the integral equation (12.121) , we rewrite 
it in the form 

Gt^siMt - s))P{s)ds = ^ {ikf (yf^P{t - s)e-"(^'-^^^«)dsj dk. 

Since (3 G -Lj^^(M-|_) with the assumed limit in (j2.16p . we apply Lebesgue's Dominated Convergence 
Theorem and compute the integrals in the principal value sense: 



/3oo r k ^, -/3oo r k^dk -/3. 



lira G'Umt-sms)ds = ^ -——-dk- . , e ^ t/2 " ,1/3 • 

The first term on the left-hand side of the integral equation (j2.12p is more tricky. First, we 
rewrite it in the form, 

Gt-s{V^{t - s))-f{s)ds = lit - s)e-^('='-*'=^o)(is) dk. 

However, if 7 S LJ^j,(M+) with the limit in (|2.16p . application of Lebesgue's Dominated Conver- 
gence Theorem yields the integral in k with a simple pole at A; = 0: 



lim fGt-s{Vo{t-s)Ms)ds = ^ [ 



dk 



_^k{k^-iVoy 

The integral is no longer understood in the principal value sense. Instead, we return back to the 
treatment of the inverse Laplace transform in (j2.5p with Re(c) > 0, use transformation o = ik, 
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and shift the contour of integration in k below the pole at A: = 0. As a result, computations are 
completed with the half-residue term at the simple pole and the principal value integral: 



hm f Gt-siVoit - s)h{s)ds = l^(-^ + r 



k'^dk \ 27, 



k^ + V^J 31^0 



Combining all computations together, we have obtained the following linear equation on /3oo and 
7oo from the integral equation ()2.12p : 

m m^i^ ~ 2|yop/3- ^ • ^ 

To deal with the second integral equation (j2.14p . we use the Fourier transform again to write 

Y r+00 ^-t{k^~ikVQ) 

and 



Ht{Vot) = — / dk 



where the integrals are understood in the principal value sense. If /?, 7 € with the limits 
()2.16p . the Lebesgue's Dominated Convergence Theorem implies that 

HtiVot), [ Ht-s{Vo{t-s))p'{s)ds^O as t ^ 00. 
Jo 

Similar to the previous computations, we prove that 

POO 

/ G'iiVot - y)uo{y)dy as t 00, 
Jo 

}}^^jyismt-s))ds = ^, 



hm / G[_,mt-s)Ms)ds 



and 



hm fGt-s{Vo{t-s))(3{s)ds = ^[ 

t~^cojQ 2.TT J_ 



3|yop/3' 

dk /3o 



^^k(k^-iV^) eiFol' 

where the last integral is computed in the principal value sense because equations (12.130 and 
()2.14p are derived in the principal value sense. 

Combining all computations together, we have obtained the following linear equation on /3oo 
and 7oo from the integral equation (j2.14p : 

- = (2-18) 
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Solving the linear system ()2.17p and (j2.18p . we obtain ()2.15p and the theorem is proved. □ 
Coming back to the original question, if ■uo(O) = 0, Mo(0) < 0, and Mq = — ^, then there is a 
finite value of to £ (0,oo) such that Ux\x=q > for all t > to. Therefore, like in the case Vq = 0, 
the function h{x, t) loses monotonicity at x = in a finite time to (where u = hx) with the only 
difference that Ux\x=o remains finite and positive as t — t- oo. We conclude that the presence of 
the advection term with Vb < in the fourth-order diffusion equation (II. 6p does not prevent the 
loss of monotonicity in x but still stabilizes the solution globally as t — )• oo. In both cases Vq = 
and Vq < 0, the monotonicity of /i in x is lost because of the inhomogeneous boundary conditions 
h \ — -1 

'''xx\x=(} — 2 ■ 



3 Solution of the original problem 



We shall now use Laplace transform ()2.ip to obtain the implicit solution to the advection-diffusion 
equation ()1.4p with a variable speed V{t). Let us denote 

W{t) = [ V{s)ds 
Jo 

and obtain the Laplace transform solution in the form: 
U{p,t) = Uo{p)e-'P'+^^'^P 

Compared with the solution (j2.4p . we have set 'y{t) = because of the third boundary condition 
in (jl.Sp . Using the inverse Laplace transform in x and recalling the definition of the Green's 
function Gt{x) (see Appendix A), we obtain the analogue of the implicit solution (|2.6p : 

u{x,t) = / Gt{x + Wit)-y)uo{y)dy-- G't_,{x + W{t)-W{s))ds 
Jo ^ Jo 

+ [ GU{x + W{t)-W{s))/3{s)ds. (3.2) 
Jo 

Now we have two unknowns /3 and W and we can set up two integral equations at the boundary 
conditions u{0,t) = and Ux{0,t) = /3{t). 

From the boundary condition u{0, t) = 0, we obtain the integral equation: 

pt roo 

- G';_s{W{t)-W{s))mds = / Gt{W{t)-y)uo{y)dy 
Jo Jo 



G't_,{W{t)-W[s))ds. (3.3) 
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To find another integral equation from the boundary condition Ux{0, t) = P{t), we differentiate 
the solution ()3.2p in x: 

roc -I pt 

/ G'tix + W{t) - y)uoiy)dy - - / G'^^ix + W{t) - Wis))ds 
Jo ^ Jo 



Ux{x,i) 

- mHt{x + W{t)) - I Ht-s{x + W{t)-W{s))l3'{s)ds 
^ Jo 

+V{t) f Gt-s{x + W{t)-W{s))p{s)ds. (3.4) 
Jo 

From the boundary condition Ux{^-,t) = we obtain another integral equation: 

/?(t)+2/3(0)i/t(VF(t)) + 2 / Ht^s{W{t)-W{s))l3'{s)ds 

Jo 

-2V{t) [ Gt-s{W{t)-W{s))l3{s)ds 
Jo 

G[{W{t) - y)uoiy)dy - [ - W{s))ds. (3.5) 



We shall prove that the system of two integral equations (j3.3p and (j3.5p determines uniquely 
the function /3(t) and V{t) locally for t > 0. The following theorem gives the result. 

Theorem 3 Assume that uq G C°°(M_|_) such that 

no(0) = 0, n^'(0) = -^, n'o"(0) = 0. (3.6) 

Then, there exists a formal solution (V, /3) of the system of two integral equations ^3. 3)) and ^3.5^ 
in the form of the fractional power series: 

oo oo 

(3{t) = /3o + 5^ /3n/4i"/', V{t) = Vo + Y. ^n/4t"/', (3-7) 

n=4 n=l 

where /3o = 'Uo(O), Vq = tiQ^^(0)/uo(0), and {/3n/4; ^{n-3)/4}J?L4 CLf^ uniquely determined. 

Proof. We substitute the series representations ()3.7p to each term of the integral equations 
(pn and (133]). It fohows from ([277]) that 

a. = ^ Vis)ds = Vot^- + ± ;r^K/4t(-^)/^ 



and 

^(„+4)/4 _ (t _ T-)("+4)/4 



If* °° 4 

6,. = ;^ l ^ Vis)ds = For^/^ + 
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Using the representation (|A.5P of the Green function with g £ C°°(M), we obtain for the three 
terms of the integral equation (13. 3p : 



f - y)uo{y)dy = g{z - at)u^{t^/^z)dz = V -4")(0)t"/^ / g(z - at)z"dz 

Jo ^^i^! ^ 

I — n ■ ■^0 JO 



k=0 " " n=2 

and 



Jo T^'^ 



At the first powers we obtain a system of hnear algebraic equations on the coefficients 

of the power series (I3.7p : 



t^/^ : -4/3o/(0) = n'o(O) / 

JO 

1 f °° 

: = -u'm g{z)z'dz, 

1 f °° 

: = -n[,"(0) / <7(^)^'d2, 



: = ^ul^\0) / - «o(0)Fo / g'{z)zdz, 







t'/' : -/34/45"(0) ^^-^^dx = l4^)(0) r g{z)z^dz - ^^(0)^0 H g'{z)z^dz 



5! - Jo 

OO 



--^[,(0)^1/4 1^ g'{z)zdz - -/(0)yo, 



and so on. Using the explicit values for the integrals ()A.9p - ()A.13P and the initial conditions p.6p . 
we obtain /Jq = ""0(0), Vq = Uq'*'' (0) /uq (0) , and the linear equation 

n'MVi/A + 8/(0) (^/34/4 + 4'^(0) + ^Fo) = (3.8) 
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Similarly, we work with the terms of the second integral equation (j3.5p : 



/■OO /"OO 1 /"OO 

/ - y)uo{y)dy = g{z - at)u',{t'/'z)dz = V u'^+'\o)t^/' / ^(^ - at)z^ 

Jo Jo ^^o""' -^0 

^o(O) E rT(-«*)' / 9^'^(^^dz + ^ -4"+i)(0)t"/4 / g{z - at)z^dz, 

k=0 ' •'^ n=l 



^0 r3/4 



I — n ■ Jo 



rat °° 1 

= / 9{z)dz = ^^5W(0)a^i 



and 



•^0 n=4 A;=0 ~'~ ''' -^0 



At the first powers we obtain a system of linear algebraic equations on the coefficients 

of the power series (|3.7|) : 



t'/' : /3o = 2^/^(0) / g{z)dz, 

Jo 

/•oo 

iV4: = 2<(0)/ <7Wzdz-4/(0), 

JO 

t2/4 . = n'o"(0) / giz)z'^dz, 







8 1 /"OO poc 







8 8 1 t \ 

t^/" : /34/4 + 7/3off(0)^i/4 - :t/3o5(0)Fi/4 = -^^O) / giz)zUz - 2u'^iO)Vo / 
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and so on. Again, using the explicit values for the integrals ()A.9p - ()A.13P and the initial conditions 
(13. 6p . we obtain /3o = u'o{0), Vq = / Uq{0) , and the linear equation 

- lu'o{0)g{0)V,/, + (^/34/4 + (0) + ^^o) = (3.9) 

The system of linear equations ()3.8p and (|3.9p has a unique solution 

^1/4 = 0, /34/4 = -ui'\0) - ^Vo, (3.10) 



/I, 



provided that 

which is confirmed. Note that the constraint Vq = u^\o)/uq{0) also follows from the contact 
equation (|1.2p obtained for sufficiently smooth solutions. Similarly, the second equation (j3.10p 
follows from the advection-diffusion equation ()1.4p after one derivative in x is taken in the limit 
X — )• and t — )• 0. 

It remains to prove that the system of linear equations obtained from the system of integral 
equations ()3.3p and ()3.5p can be solved at each order of i"/^ for n G N. From the previous 
computations, we can deduce that the first integral equation at t^^+^^Z^ gives a linear equation 
on variables (/3„/4, V(„„3)/4) of the power series p.7p : 

- /3„/4/(0) ^ ^3/^ cirr + -_n[,(0)y(„_3)/4 g'{z)zdz = ■■■ , (3.11) 

where the dots denote the terms expressed through derivatives of uq{x) at x = and the previous 
terms of the power series (j3.7p . Similarly, the second integral equation at t"/^ gives another linear 
equation on variables (/3n/4) ^(n-3)/4)- 

/3„/4 - ^/3o5(0)V(„_3)/4 = • • • . (3.12) 
The system of linear equations (13. lip and (I3.12p is non-degenerate if 

The left-hand side {C„}„gN is computed numerically (see Figure[2]). It is a monotonically increas- 
ing sequence that approaches closely to 1 at rz = 8, where Cs ~ 0.96, and n = 9, where Cg « 1.04. 
Therefore, the linear system is non-degenerate and a unique solution for {f3n/4,V(^n-3)/4) exists 
for any n G N. □ 
In the present time, we cannot prove yet that the system of integral equations (j3.3p and 
(|3.5p leads to a finite-time blow-up, according to the conjecture in [1]. Nevertheless, numerical 
computations show that the blow-up holds for a generic set of initial data. Figure [3] shows the 
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Figure 2: Numerical approximations of Cn defined by (|3.13p . 



behavior of functions (3{t) and V{t) near the blow-up time. It follows from this figure that 
f5{t) = hxx\x=o — ?• at the same time as V{t) — )• — oo with — )• Co, where Co > is 
a numerical constant. In other words, we conclude with the conjecture that /3(t) ~ V{t)~^/'^ as 
V{t) — >• — oo in a finite time to G (0,oo). 




1000 10000 100000 lOOOUOO lUOOOIKIO 

T (log. scale) 



Figure 3: Numerical computations of (3{t) and V{t) for the advection-diffusion equation (jl.ip . 
We thank the authors of [IJ for this numerical figure. 
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A Green's function 



Let us define the fundamental solution of the fourth-order diffusion equation: 

ht + hxxxx = 0, X £ M, t > 0, ('A Ti 

h\t=o = d{x), xGM, ^ 

where (5 is a standard Dirac delta-function in the distribution sense. The fundamental solution 
is usually referred to as Green's function and we shall denote it by 

h{x,t) = Gt{x), xGM, tGM+. 

Using the Fourier transform in x, we can obtain the explicit expression for Green's function: 

Gt{x) = — / e~"''+'^''dk = - / e'"'\osikx)dk. (A.2) 



2tt 7_oo vr 
In particular, we have Gt{—x) = Gt{x) for all x G M and 







G;'(o) = -i/°°.v-*d* = -EH§, ,A.4) 

where T is the standard Gamma function. The Green's function can be represented in the self- 
similar form by 

1 / X \ 1 4 

Gtix) = [^^) , g{z) = -J^ e-^\os{kz)dk, (A.5) 

where g G L2(M) n L~(M). Therefore, Gt decays to zero as t — )• oo in any norm for p > 2. In 
particular, < ||5'||L°°/i"^'''^; |G'f(x)| < \\g'\\L°° /t^^"^ , and so on, for any x G M. 

By the stationary phase method (see, e.g.. Chapter 5 in [3]), g{z) and all derivatives of g{z) 
decay to zero as |z| — >• oo faster than any algebraic powers. This gives the decay of Gt{x) and any 
x-derivative of Gt{x) as |x| — )• oo for any fixed t > 0. Although Gt and g are not L} functions, 
they satisfy the normalization conditions: 

Gt{x)dx = I g{z)dz = 1, t > 0. (A.6) 

The even function 5^ : R — t- M satisfies the ordinary differential equation 

subject to the initial values 
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and the decay behavior as |2;| — oo. It is clear from the differential equation that g G C°°(M) 
satisfies a number of integral constraints: 
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